MATHEMATICAL ENGINEERING TECHNICAL REPORTS Hamiltonian Monte Carlo with Explicit, Reversible, and Volume-preserving Adaptive Step Size Control
نویسندگان
چکیده
Hamiltonian Monte Carlo is a Markov chain Monte Carlo method that uses Hamiltonian dynamics to efficiently produce distant samples. It employs geometric numerical integration to simulate Hamiltonian dynamics, which is a key of its high performance. We present a Hamiltonian Monte Carlo method with adaptive step size control to further enhance the efficiency. We propose a new explicit, reversible, and volume-preserving integration method to adaptively set the step sizes, which does not violate the detailed balance condition or require a large increase in computational time.
منابع مشابه
Adaptive Step Size for Hybrid Monte Carlo Algorithm
We implement an adaptive step size method for the Hybrid Monte Carlo algorithm. The adaptive step size is given by solving a symmetric error equation. An integrator with such an adaptive step size is reversible. Although we observe appreciable variations of the step size, the overhead of the method exceeds its benefits. We propose an explanation for this phenomenon.
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Graver Basis for an Undirected Graph and Its Application to Testing the Beta Model of Random Graphs
In this paper we give an explicit and algorithmic description of Graver basis for the toric ideal associated with a simple undirected graph and apply the basis for testing the beta model of random graphs by Markov chain Monte Carlo method.
متن کاملReliability and Sensitivity Analysis of Structures Using Adaptive Neuro-Fuzzy Systems
In this study, an efficient method based on Monte Carlo simulation, utilized with Adaptive Neuro-Fuzzy Inference System (ANFIS) is introduced for reliability analysis of structures. Monte Carlo Simulation is capable of solving a broad range of reliability problems. However, the amount of computational efforts that may involve is a draw back of such methods. ANFIS is capable of approximating str...
متن کاملA New Approach for Monte Carlo Simulation of RAFT Polymerization
In this work, based on experimental observations and exact theoretical predictions, the kinetic scheme of RAFT polymerization is extended to a wider range of reactions such as irreversible intermediate radical terminations and reversible transfer reactions. The reactions which have been labeled as kinetic scheme are the more probable existing reactions as the theoretical point of view. The ...
متن کاملGeometric integration of conservative polynomial ODEs ✩
We consider systems whose Hamiltonian is of the form H(q,p) = 2p + V (q), where the potential V is either cubic or quartic with no cubic terms. For most of these systems (in the measure sense) we give an explicit numerical integration method that preserves both phase space volume and the value of the Hamiltonian. This is exemplified in the Hénon–Heiles system. An application is to the hybrid Mo...
متن کامل